General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



N77-15809 


(NASA-TM-X-73196) COMMENT ON THE RELATION 
BETWEEN THE NONADIABATIC COUPLING AND THE 
COMPLEX INTERSECTION OF POTENTIAL ENERGY 
CURVES (NASA) 20 p HC A02/MF A01 CSCL 20H Unclas 

G3/72 11537 


NASA TECHNICAL 

MEMORANDUM 


NASA TM X-73,196 



COMMENT ON THE RELATION BETWEEN THE NONADIABATIC 
COUPLING AND THE COMPLEX INTERSECTION OF POTENTIAL 
ENERGY CURVES 

Ricliard L. Jafte 

Ames Research Center 
Moffett Field, CaUf. 94035 




1. Rtport No. 2. GoMrnitMnt Accwsion No. 

NASA TM X-73,196 


4. Title and Subtitle 

COMMENT ON THE RELATION BETWEEN THE NONADIABATIC 
COUPLING AND THE COMPLEX INTERSECTION OF POTENTIAL 
ENERGY CURVES 


7, Author{s) 

Richard L. Jaffe* 


9. Performing Orgeniutlon Name and Address 

Ames Research Center, NASA 
Moffett Field, Calif. 94035 


t2. Sponsoring Agency Name and Addreu 

National Aeronautics and Space Administxation 
Washington, D.C. 20546 


15. Supplementary Notes 

*NASA/NRC Research Associate* 1974-1976. 


3. Recipient's Catalog No. 


5. Report Date 


6. Performing Organization Code 


8. Performing Organization Report No, 

A-6705 


to. Work Unit No. 

506-16-12 


n, Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Memorandum 


14. Sponsoring Agency Code 


16. Abstract 


Simple relations are discussed that provide a correspondence between 
the complex intersection of two potential surfaces and the nonadiabatic 
coupling matrix element between those surfaces. These are key quantities 
in semiclassical and quantum mechanical theories of collision- induced 
electronic transitions. Within the two-state approximation, the complex 
intersection is shown to be directly related to the location and magnitude 
of the peak in the nonadiabatic coupling. Two cases have been considered; 
(I) the avoided crossing between two potential surfaces, and (II) the spin- 
orbit interaction due to a halogen atom. Comparisons are made between 

the results of the two-state model and the results of ah quantum 

chemical calculations. 


17. Kay Words (Sugg«tad by Auttwr(sM 

Atomic and molecular collisions 
Electronic transitions 


19. Sacurity Qa«lf. (of this raporti 

Unclassified 


18. Distribution Statament 


Unlimited 


STAR Category - 


20. Sacurity Ctassif. (of this paga) 

21. No. of Pages 

Unclassified 

19 



*For sala by tha National Tachnical Information Sarvica, Springfiald, Virginia 22161 





















COMMENT ON THE RELATION BETWEEN THE NONADIABATIC COUPLING 


AND THE COMPLEX INTERSECTION OF POTENTIAL ENERGY CURVES 

Richard L. Jaffa 
Ames Research Center 

SUMMARY 

Simple relations are discussed that provide a correspondence 
between the complex Intersection of two potential surfaces and the 
nonadlabatlc coupling matrix element between those surfaces. These 
are key quantities in semiclasslcal and quantum mechanical theories 
of collision-induced electronic transitions. Within the two-state 
approximation, the complex intersection is shown to be directly 
related to the location and magnitude of the peak in the nonadlabatlc 
coupling. Two cases have been considered: (I) the avoided crossing 

between two potential surfaces, and (II) the spin-orbit Interaction 
due to a ^P halogen atom. Comparisons are made between the results 
of the two-state model and the results of db ■initio quantum chemical 
calculations. 


1. Introduction 


The first quantum mechanical calculations of collision- induced 
electronic transitions for model trlatomlc systems have been reported recently 
[1,2]. These studies considered electronic-to-vibrational [1] and 
electronic-to-rotational [2] energy transfer in nonreactive collisions of 
halogen atoms with hydrogen molecules. In addition, semiclassical 
calculations have been carried out for electronic transitions in h"^ 
and H 2 +F(^Pi/ 2 ) [^3 collisions based on the method developed by Miller 
and George [5]. Previously, quantum mechanical [6] and semiclassical [7] 
calculations had been reported for fine-structure transitions in atom-atom 
collisions. In the quantum mechanical studies , electronic transitions are 
induced by the nonadiabatic coupling matrix elements between the initial 
and final adiabatic electronic states (potential coupling) . These coupling 
matrix elements can be calculated from the electronic wave functions or , to 
a good approximation, from the adiabatic electronic potential surfaces them- 
selves [1,8]. For diatomic systems, adiabatic potential surfaces of the 
same symmetry and multiplicity do not intersect for real values of the 
nuclear coordinates, but complex-valued intersection points (branch points) 
do exist [9-13]. For larger systetis, the surface of complex intersection 
always has a higher dimensionality than the surface of real intersection. 

The semiclassical theory makes use of the complex intersections to effect 
transitions between the potential surfaces. The collisions are simulated 
by classical trajectories which are analytically continued from one potential 
surface to the complex intersection and back to the second potential surface. 
The probability for making such a transition shows a roughly Inverse dependence 
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on the distance from the real axis to the complex Intersection point. The 
two methods are similar in that the factors which govern the magnitude of 
the transition probability are derived from the adiabatic potential surfaces 
involved. Davis and Pechukas [14] and Pechukas et al. [15] have shown that 
the semiclassical treatment is indeed a limiting case of the quantum mechanical 
theory. However, the problem of locating the complex Intersection points 
must still be solved before semiclassical calculations of electronic transitions 
can be undertaken. 

Ah initio molecular orbital calculations of potential energies at 
complex values of the nuclear coordinates have recently been performed 
for four systems. The location of the complex intersection points has been 
determined for the 3a and 4a states of HeH [10], the two lowest singlet 
states of [11] and the lowest and states of collinear H^F [12] 
and HF"* [13]. In the first two systems, the potential surfaces for the two 
states exhibit avoided crossings for real-valued geometries, while, in the 
latter two systems, the potentials are parallel in the limit of separated 
+ F and + F, and are split by the spinr-orbit coupling of the fluorine 
atom. There is no avoided crossing of the HF*^ potential curves which diverge 
as the internuclear separation is decreased. For the collinear H^^F system, the 
two states show (in addition to the asymptotic behavior discussed above) 
an avoided crossing at small HF distances. These studies have shown qualita- 
tive relationships for the location of the avoided crossing, the actual 
complex Intersection point, and the extremvnn in the nonadlabatic coupling. 

For these systems, the complex intersection points were located by the analytical 
continuation of rational- fractions which had been numerically fit to the 
energy difference between the potential surfaces [16,17] . The real part of 


the complex nuclear coordinate at an intersection point was found to be 
approximately equal to the value of that coordinate at both the extremum 
in the nonadlabatlc coupling and the avoided crossing. 

In this paper, the relationship between the complex intersection point 
and the nonadlabatlc coupling is explored. It is shown that formal 
connections exist for special cases within the two-state approximation. 

The results for these cases are presented in Sections 2 and 3 and discussed 
in terms of the results of the complex-valued db initio molecular orbital 
calculations . The above relationship can be used to help locate the complex 
intersection points for cases when the results of complex-valued molecular 
orbital calculations are not available. Section 4 contains some concluding 
remarks. 

2. Case I: The avoided crossing 


In the two-state model of an avoided crossing, the electronic Hamiltonian 
matrix is given by 


H 


Hii(r) Hi2(r) 


H (r) H„„(r) 

12 ~ 22 ~ 


( 1 ) 


where r is a vector that represents all the nuclear coordinates. The 
potential energies of states 1 and 2 are 


Ej - (-1)^ f , 


where 


AE 


[(H 22 " > 4Hf^3 


1/2 


( 2 ) 


(3) 
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is the energy difference between the two adiabatic states. The nonadiabatic 
coupling between the states for coordinate a, In this approximation » Is 
given by [8,12] 


»12 - » 11 > 
aa 



In the following treatment, we assume that the electronic transition 
occur ^ only through motion along one generalized coordinate r and that 
the coupling matrix element is Independent of that coordinate. We also 

define the quantity f equal to and use primes to signify partial 

differentiation with respect to r. 

Extrema in AE occur at those values of r for which 

AE' » ff'/AE , (5) 

that is, when f = 0 or f * 0. Examination of AE” reveals that AE 
is always a minimum when f - 0. This case corresponds to the avoided inter- 
section ( r =» r ) . ; Extrema in the nonadiabatic coupling are found by solving 


H,„f” 

= — 0 . ( 6 ) 

AE^ AE** 

c 

From eq. (6) we conclude that the system will exhibit an extremum in the 
nonadiabatic coupling at a point of avoided Intersection only if f = f” = 0. 
To determine the condition for r^^^ = Rer^, the real part of a point of 
complex intersection, we first expand f (r) in a Taylor series about r^; 



(r - rj2f«(rj (r - r J^f»”(r ) 
f(r) = (r - r^)f-(r,) + . . . 


( 7 ) 
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At a complex intersection point (AE ■ 0) fCr^) must be pure imaginary and 

equal to ±i2H^2* Since f is a real function of the complex variable r, 

Re Tq ■ r^jj only if all even power terms in eq. (7) are zero (i.e., f^^^^(r^) = 0) 

This condition encompasses the one found above for r = r . 

c m 

Two functional forms for f which exhibit the property Rer^ “ 


are 


and 


f = ar - b 


f = A sin(ar - b) 


( 8 ) 


(9) 


a, b, and A can be functions of any of the coordinates except r. Equation (8) 
has been used in the Landau-Zener treatment of electronic transitions [18] . 

In both cases, 

r = b/a and AE(r ) = 2H,„ . 

m m i 2 

The location of the complex intersection points for f given by eq. (8) 


are 


iAE(r^) 

m 

r ± 

m a 


( 10 ) 


and the magnitude of the nonadiabatic coupling at r^ is a/2AE(r^) . Thus 


|Xj 2 I ^ single maximum of height 1/2 Imr^. This provides a demonstra- 

tion of the behavior found in the semiclassical calculations; namely, the 
transition probability falls off inversely with the distance from the real- 
axis to the complex intersection point. When f is given by eq. (9) , 
complex Intersection points are given by 


r ± slnh • I — ■: 1 

m a L ^ J 


( 11 ) 
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and 




Aa 

2AE(r ) 
' m 


( 12 ) 


while !x^2r^ simply proportional to 2/lmrQ, it still decreases with 

increasing Imr^. for both cases, the strength of the transition is proportional 

to the slope of f at r - r . 

x& 

The half width, W, of the peak in Xj 2 ^ Pi^o^ldes a measure of the 
degree to which the electronic transition is localized around r^. W is 
equal to r - r , where r“ are obtained by inversion of eq. (A) with 

Xi2r “ ^i2r^^c^^^' ^ 

r” = r ± AE(r )/a , (13) 

m m 

so the half width of the peak in the nonadiabatlc coupling matrix element 
is simply 2AE(r^)/a * stronger the coupling, the more 

localized the transition. ^ 

The reliability of this model can be tested by comparison with the 
db 'initio data for HeH^ and For HeH^ [10] , the energies of the 3a and 4a 

states (E„ and E, ) were determined from the third and fourth eigenvalues of a 

3a 40 

16 X 16 electronic Hamiltonian matrix as a function of r, the internuclear 
separation. The avoided crossing was located at r^ = 3.64A0 bohr, the value 

-f* 

of r that minimized AE « E ' - E„ . The electronic transition in H. + H + H 

collisions occurs between the first and second singlet potential surfaces of 
H^. In this system, for colllnear or Isosceles triangle geometries, there are 
two nuclear coordinates,* r and R, and the loci of avoided crossings and complex 


* For colllnear geometries, r and R are the distances between the central atom 
and the two outer atoms with r < R. For isosceles triangle geometries, 
r and R are the base and altitude of the triangle, respectively. In the 
first and second singlet states of r corresponds to the vibrational 
coordinate for and respectively, and R corresponds to the translational 
coordinate. 
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• I * 1 r j I / 

Intersection points lie on curves of (r,R) through configuration, space. 

Ab irvitio calculations were performed to locate r^ and r^ for fixed values 

of R [11], It was found that r^ is roughly unchanged for values of 

R > 6 bohr, indicating that the electronic transition is effected through 

vibrational motion during the collision. ’ 

For these two systems, a can be obtained from the slope of a linear f:t 

of r versus f = ±{[AE(r)]^ - [AE(r )]^}^/^** and Imr is predictea from 

m o 

eq. (10). The results are shown in table 1 for four sets of data from the 

ah initio calculations. For small errors in the results of the model 

are less than 1 / 2 % when compared with the dh initio data, as can be seen for 
|-[ ^ 

the HeH and H 3 R = 9 bohr results in table 1. However, for larger di 2 > 
the model breaks down. f(r) is still nearly linear and Imrp can be predicted 
to within 90% of the true value, but the model cannot account for the shift in 
Rer^ away from r^. The fitting of f(r) to a rational-fraction in r [17] 
improves the two-state model. However, r^ and r^ could no longer be determined 
analytically and the relation ^jjj “ ® r^ would no longer hold. 

Thus, this simplified model of an avoided crossing is quite accurate when 
applied to systems with weak coupling ( 11 ^ 2 ) which a definite relation 
exists between the factors that influence the quantum and semlclasslcal 
transition probabilities . This implies that the electronic transition is 
localized at r ~ add the complex intersection points can be located with 
the use of the real-valued potential energy surfaces. However, for stronger 
coupling, the correspondence between the two methods no longer holds, and 


** f is positive for r < rjj^ and negative for r > r^. 
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one must know the potential energies at complex values of the nuclear 
coordinates to carry out semlclasslcal calculations of electronic transitions.. 

3. Case II: Spinr-orbit coupling and fine-structure transitions 

The splnr'urbit interaction for atom-atom and collinear atom-^iatora collisions 
involving halogen atoms can often be treated within the framework of a 
two-state approximation [7,12,19,20], In these systems, the nonrelativ.istic 
description of the halogen atom has a degenerate ground state which splits 
into ^^ 1/2 upon the inclusion of spinr-orbit coupl5.ng. 

The potential curves for an A + X collision (A is any nondegenerate atonic 
or diatomic species and X is a halogen atom) are parallel at large A-X 
separations with an energy difference equal to X, the halogen atom spin orbit 
splitting. These asymptotes correspond to A + X (^^ 3 / 2 ) A-X 
The potential curves are subject to the same constraints regarding real and 
complex intersections as in the case of an avoided crossing. 

Within this approximation, the complete Hamiltonian matrix is given by 
eq. (1) with H, , = + X/3,H^^ “ V„, and » (2)l/^(X/3). and V are 

the nonrelativistic potential energy surfaces for the lowest H and E states 
of the system [12]. For this case, eq. (3) becomes 

AE « [(AV + X/3)^ + 8X^/9]^^^ , (14) 

where AV = V_ - V . At the asymptotic limit, AV ® 0 and AE equals X. 

li u 

The complex intersections occur when AV * AVq , 

AVq * X exp [±i cos ^(-1/3)] • (lb) 
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If AV + (X/3) equals zero, there will still be an avoided crossing with its 
associated complex intersection points as in case I. However, additional 
complex intersection points are caused by the asymptotic behavior of AV. 
Since AV can be qualitatively represented by functions like x“’^ or e“^, 
the relations found for case I are not applicable to the spin-orbit case, 

Preston, Sloane, and Miller [7] have fit AV to an exponential function 
in the internuclear separation, r, 

AV(r) = A exp (-otr) (16) 


for the Xe + F system in which V^ is the ground state. We use this form 
in the following treatment of fine-structure transitions. For collinear 
+ F, r represents the HF separation with A and a taken to be functions 
of the HH bond length. Using eqs. (15) and (16), the complex intersection 



Thus, Im r^ depends only on the difference in shape between the nonrelativistic 
potential curves while Re r^ depends on both the strength of the spin-tjrbit 
coupling and the nature of AV, The faster AV goes to zero with increasing 
r, the greater the electronic transition probability, since the complex 
intersection is located closer to the real axis. 
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Tha nonadlabatlc coupling, as given by eq. (4), has extrema at the 
values of r for which 



These results, which are summarized in table 2, are also similar to the 
relations found for case I. However, Im r^, ^ depend 

on the strength of the spin-orbit coupling. 

The available db initio data can be used to test this model as was done 
for case I. Calculations of and as a function of r, the H-F 
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separation t have been performed for H^F [12] in collinear geometries with the 
H-H bond fixed at 1.4 bohr (the equilibrium separation in H^) and HF"^ [13,20]. 
These calculations incorporated the spin-orbit interaction within the framework 
of the two-state model discussed above. In the former system, is the 
ground-state surface at small r, but a shallow well in the H-state causes 
V„ < V„ for r greater than 4.7 bohr. Thus, AV exhibits a minimum of 
-2.8 X 10 ** hartrees at 5.4 bohr and cannot be represented by eq. (16). 

However, r^ can be predicted from r^ and Xj^2r^^c^ using eqs. (17) and 
(20). The resulting complex intersection point, 2.023 ± 0.6331 bohr, is in 
fair agreement with the db inttio value, r^ = 2.041 ± 0.4291. For HF , 

V < V for all r, and AV is better represented by eq. (16). However, a 

ns 

plot of In AV versus r shows considerable curvature [13,20]. The data from 

the db initio MO calculations of AV are in good agreement for r < 5 , but 

the results of Jaffe et al. [13] fall off less steeply than do those of ref. [20] 

at larger r. The calculated values of r^ and X]^2r^’^c^ various fits 

of these data to eq. (16) are presented in table 3. It can be seen that good 

values for Im r and x, (r ) are obtained from the fit to the r >6 
0 2 2ir c 

data and that no fit reproduces Re end r^ accurately. However, use of 
the ab initio data for r^ and leads to an accurate prediction of the 

location of the complex intersection point, even though the aobiuii data points 
oannot be aoQuratety fit to eq. (16), 

For HF'*’, it has been noted [7] that AV could be better represented by 

AV = A exp[-ar + 8r^] . (22) 

In this case, eq. (18) becomes a transcendental equation in r and an 
analytic solution for r^ is not possible. Numerical determinations of the 

. . 12 ... 


extrema of have been carried out for various values of A, ot, and 3. 

The location of the complex Intersections have been shown to be [7,13] 

to - (“" - «[ln (f\ ± 1 cor* (+ • C23) 

The relations between r^ and very close agreement with the 

above case for 3 < 0.03 a, but only at 3 = 0 is r^ exactly equal to Re r^. 
In addition, an empirical relation has been found between ”0 


(24) 

In ref. [13], it is reported that is quite sensitive to which data 

points are used to determine A, a, and 3 by a least-squares fitting procedure. 
The empirical relation given by eq. (21) suggests that the best value of 
Im tjj is that obtained from the parameters for which the fitted and ab tnitio 
data for Xj^2r^^c^ are in closest agreement, provided, of course, that 
3 < a^/4 In (A/X). This criterion provides a means of selecting the optimum 
parameters in eq. (22) when an ab initio calculation of the complex intersection 
is not available. Examination of the detailed results of the least-squares 
fits used in ref. [13] confirms the validity of this procedure. However, r^ 
can only be located to within (0.03 + O.Oli) bohr by this method. 

4. Conclusions 

It has been shown that simple relations exist among the quantities which 
govern the electronic transition probability in quantum and semiclassical 
scattering theory for two special cases. For an avoided crossing, these 

: 13 ^ 




relations hold when the coupling between the electronic states small 

and the transition Is localized at the avoided Intersection. This condition 
should encompass most cases with large transition probabilities. In seml- 
classlcal calculations > the complex intersection points can be located with the 
use of the slope of a linear fit of this wayi semiclassical 

calculations can be carried out using only real-valued potential energy 
surfaces as Input. 

For the spln-orblt case, it has been shown that the nonadiabatlc coupling 
and location of the complex Intersection points are very sensitive to 6V, 
the difference between the nonrelatlvlstlc potential surfaces. Simple relations 
exist for these qviantitles in systems with relatively short-range interactions 
(large a's and small or negative fi's in eqs. (16) and (22)). 

Again, this provides a method to help locate the complex intersection 
points needed for semiclassical calculations of electronic transitions when 
only real-valued potential energy surfaces are available. However, in this 
case more caution is needed as the location of complex Intersection points 
is very sensitive to the functional form used to fit the AV versus r data 
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Comparison of predictions from the two-state model of avoided crossings 
with complex-valued ab initio MO data 


Case 

r 

m 

(bohr) 

AE(r ) 
(hartree) 

a“^ 

(hartree 
bohr ^) 

^0 

(bohr) 


r^ predicted^^ 
(bohr) 

quality*^^ 
of fit 

^.d) 

HeH 

3.6440 

0.00294 

-0.0784 

3.6440 ± 

0.03751 

3.6446 ± 0.03751 

1.48 

H** collinear®^ 

3 

2.447 

0.00117 

-0.0542 

2.4482 ± 

0.02181 

2.449 ± 0.02181 

1.92 

(R - 9 bohr) 








collinear®^ 

2.448 

0.01683 

-0.0527 

2.3926 ± 

0.29591 

2.446 ± 0.31931 

5.39 

(R = 6 bohr) 








+ e) 

H3 triangular 

2.468 

0.02621 

-0.0515 

2.3464 ± 

0. 44391 

2.468 ± 0.50901 

4.67 

(R = 6 bohr) 









a) Determined by linear least squares fit of f(r). 


b) r„ predicted = r ± iAE(r_,)/a. 
0 m m 


c) Defined as 100 x 


± i 




d) Data from ref. [10], 

e) Data from ref. [11] for [5s, 2p]-^ (3s, Ip) basis set. 


Table 2 

Summary of Interrelations for case 11 

W - 2.98199/a. 

W - 2.19722/a, V^. > 

IXi2r<’'c)i - 0-35355 /c. V,. > 

IW'c’l ■ ''n" '';. 

Im r -±1.91063/a, V„ > V„ 

0 ' “ ^ 

Im r^j - ±1.23096/a, V^. > Vj^ 


Table 3 

Conq)arison of the model ab initio MO data for HF 


Case*^ 

a 

(bohr 

A 

(hartree) 

r 

c 

(bohr) 

Xi2r(^) 

(bohr"'^ 

r 

0 

(bohr) 

fit to r < 5 data 

0.91 

o.7e\v^ 

6.619 

0.322 

6.619 ± 1.353i 

fit to r > 6 data 

0.39 

0.043 

"^^05 

0.138 

8.305 ± 3.1561i 

fit to all data 

0.62 

0.235 

7.822\^ 

0.219 

7.822 ± 1.985i 

cd> initio results 

- 

- 

7.4 

0^141.,^^ 

7.495 ± 3.009i 

(!) predicted from 






'c* >‘l2r<V 

0.399 

0.035 

- 

- 

\7.4 ± 3.087i 

(ii) predicted from 





\ 


0.409 

0.039 

7.495 

0.145 



a) ab initio data taken from ref. [13]. 






